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Abstract 

We investigate theoretically the feasibility of generation and am- 
plification of terahertz radiation in aligned achiral carbon nanotubes 
(zigzag and armchair) in comparison with a superlattice in the pres- 
ence of a constant (dc) and high-frequency (ac) electric fields. The 
electric current density expression is derived using the semiclassical 
Boltzmann transport equation with a constant relaxation time with 
the electric field applied along the nanotube axis. Our analysis on 
the current density versus electric field characteristics demonstrates 
negative differential conductivity at high frequency as well as photon 
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assisted peaks. The characteristic peaks are about an order of mag- 
nitude styronger in the carbon nanotubcs compared to supcrlattice. 
These strong phenomena in carbon nanotubes can be used to obtain 
domainless amplification of terahertz radiation in carbon nanotubes at 
room temperature. 

1 Introduction 

Terahertz (THz) region of the electromagnetic (EM) spectrum refers to the 
frequency range between the mid infrared and the microwave region (0.3 THz 
- 10 THz). This EM region has a lot of promising applications in various 
areas of science and technologies include astronomy, broadband communi- 
cation, pollution monitoring , poison gas sensing, etc. [1, 2j. In spite of 
the potential applications, this region of the EM spectrum is yet to be fully 
exploited, due to the limited availability of device sources and detectors 
[3]. Over the past decade, terahertz science and technology has advanced 
considerably with both optical-bench-based systems and solid state quan- 
tum cascade lasers [3] . The main drawback of these quantum cascade THz 
sources is low working temperature, making it difficult to maintain popu- 
lation inversion at room temperature [4]. So far the highest temperature 
attained for a THz source operating at 3 THz frequency is 160 K [5] . Creat- 
ing a compact reliable source of THz radiation that could operate at room 
temperature still remains one of the most formidable challenges of contem- 
porary applied physics. Most recently, carbon nanotube has been identified 
for THz technology due to its superb physical properties. Different propos- 
als of carbon nanotubes for THz applications have been made; that ranges 
from frequency multipliers [U [7] , THz amplifies [8] , switches [9] to antennas 

Carbon nanotube (CNT) is single layer nanometer-size cylinder made 
of carbon atoms. The cylinder has no bulk atoms, only surface atoms. It 
is a layer of graphene seamlessly rolled into cylinder. CNTs have intrigu- 
ing physical properties (mechanical, electrical, optical, etc.) due to their 
unique structure; combine with of one-dimensional solid state characteris- 
tics with molecular dimensions [11]. These properties vary with each CNT 
fundamental indices, (n, m) which specify the diameter and wrapping angle 
as the graphene sheet seamlessly wrapped into a CNT. As n and m vary, 
conduction ranges from metallic to semiconducting [12] , with n = m been 
metallic and m = been semiconducting. The band gap of the semicon- 
ducting CNTs which is inverse diameter (d) dependent, ranges from 0.2 eV 
to 2 eV with E g = 0.9/dnm [13], where d = a (n 2 + m 2 — mnj^/f is 
the nanotube diameter and a Q = 0.246 nm is graphene lattice constant. 
The peculiarity of CNTs (e.g., strong nonparabolic dispersion law) raises 
the question of whether these structures could be used for high-order har- 
monic generation |15j . Different theoretical models and experimental 
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techniques are being pursued to demonstrate the feasibility of CNT for THz 
applications[16j - |20j. These ranges from THz generation Cherenkov-type 
emitters base on CNTs and hot electrons in quasimetallic CNT (n — m = p 
with p as a non-zero interger), frequency multiplication in chiral-nanotube- 
base superlattices controlled by a transverse electric field, and THz radiation 
detection and emission by armchair nanotubes in a strong magnetic field [16j. 
ab initio molecular dynamic simulations [TTJ. CNT capacitor circuit model 
|18j . impedance transmission model [19] . effective conductivity model [20] . 
just to mention a few. 

Nonlinearity of the current density characteristics of CNTs gives rises 
to a generation of harmonics of microwaves and direct current (dc) gener- 
ation as well as terahertz generation [TJ [2"2"1, 123] . However, Bloch oscilla- 
tions of electrons within a CNT cause an appearance of negative differential 
conductivity (NDC) for dc fields larger than the critical electric field |15| . 
This leads to electric current instabilities and formation of domains in the 
CNTs [23]. These current instabilities are destructive to the formation of 
the terahertz frequency. Therefore, it is important to consider the scheme 
of generation of terahertz radiation than can suppress electric instabilities 
while preserving high-frequency gain at room temperature. Theoretical re- 
alizations of superlattice oscillators (SL) with ac bias have been reported 
|25] - [27], that demonstrated domainless amplification of terahertz radia- 
tion. However, there are limited reports of such effects in CNTs, especially 
simultaneous use of dc-ac electric field. We present here, a report on a 
theoretical investigation of zigzag single wall carbon nanotube (z-SWCNT) 
and armchair single wall carbon nanotube (a-SWCNT) in comparison with 
superlattice (SL) under the influence of a time dependent dc-ac electric field 
for THz frequency generation and amplification at room temperature that 
indicates an order of magnitude effective suppression of the electric insta- 
bility in the one dimensinal CNTs in comparison with the two-dimensional 
SL as reported in refs [25] - [27] . 

2 Theoretical model 

Following the approach of refs. [IB]- [13], we consider an undoped a-CNT 
(i.e., zigzag and armchair ) exposed simultaneously to a dc-ac electric field. 

E(t) = E Q + E lC osujt (1) 

where E Q is the constant electric field, E\ and oj are the amplitude and 
the frequency of the ac field respectively. The investigation is done within 
the semiclassical approximation in which conduction electrons with energy 
below the energy of the interband transitions move in the crystalline lattice 
like free quasi-particles with dispersion law extracted from quantum theory. 
Taking into account the hexagonal crystalline structure of a rolled graphene 
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in a form of CNT and using the tight binding approximation, the energy 
dispersion for zigzag and armchair CNTs are expressed as in Eqns. @ and 
Q, respectively pQ. 



£(sAp (p ,p z ) = £ s (p z 



7o 



1 + 4cos(ap z )cos(-^=sAp ip ) + Acos 2 (-^=sAp lf 



1/2 
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1 + 4:Cos(ap z )cos(— =sAp ip ) + Acos (— =sAp ip ) 
v3 v3 



(3) 



where 70 ~ 3.0eV is the overlapping integral, p z is the axial component 
of quasi-momentum, Ap^ is transverse quasi-momentum level spacing and 
s is an integer. The expression for a in Eqns. ([2]) and ^ is given as 
a = 36/27T, b = 0.142nm is the C-C bond length. The — and + signs 
correspond to the valence and conduction bands respectively. Due to the 
transverse quantization of the quasi-momentum, its transverse component 
can take n discrete values, p v = sAp^ = (7r^/3s)/an (s = l...,n). Unlike 
transverse quasi-momentum p v , the axial quasi momentum p z is assumed 
to vary continuously within the range < p z < 2n/a , which corresponds to 
the model of infinitely long of WNT (L = 00). This model is applicable to 
the case under consideration because of the restriction to the temperatures 
and/or voltages well above the level spacing [U [12], ie. ksT > £c, AS, 
where fee is Boltzmann constant, T is the temperature, £c is the charging 
energy. The energy level spacing A£ is given by 

Af=^, (4) 

where vf is the Fermi velocity and L is the carbon nanotube length [5] 

Employing Boltzmann equation with a single relaxation time approxi- 
mation 

g/fo*) { cE(t) df(P,t) = [f(p,t)-fo(p)] _ 
dt dp t 

where e is the electron charge, fo(p) is the equilibrium distribution func- 
tion, (p, t) is the distribution function, and r is the relaxation time. The 
electric field E is applied along CNs axis. In this problem the relaxation 
term r is assumed to be constant. The relaxation term of Eqn. §5§ de- 
scribes the effects of the dominant type of scattering (e.g. electron-phonon 
and electron- twistons) . For the electron scattering by twistons (thermally 
activated twist deformations of the tube lattice), r is proportional to m 
(armchair dual index (m, m)) |27| and the I — V characteristics shows that 
scattering by twistons increases and decreases in the negative differential 
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conductivity (NDC) region; the smaller m is, the stronger this effect. Quan- 
titative changes of the I—V curves turn out to be insignificant in comparison 
with the case of r = const [27J EE] • The distribution functions for zigzag 
CNT could be expanded in Fourier series [1] as, 

n 

f {p) = A Pip Y, S(p v - aA Pv ) frse ibrp % (6) 

n 

fo(p, t) = A Pip ]T S(p v - sA Pv ) £ f rs e lbr P^ v (t), (7) 

5=1 r^O 

Here, the energy unit h = 1 has been used, 5(x) is the Dirac delta func- 
tion, r is summation over the stark component, f rs is the coefficient of the 
Fourier series and & v is the factor by which the Fourier transform of the 
non-equilibrium distribution function differs from its equilibrium distribu- 
tion counterpart. All other parameters are as defined above. f rs is expressed 
as. 

frs = ~ !^ d Pz (8) 

2-n-AppS J 1 + exp(£ s (p z )/k B T) 

Substituting Eqns. ^ and ([7]) into Eqn. ([5]), and solving with Eqn. ([!]) we 
obtain 

fc=oo m=oo v ' 

where [3 = eaE/h, Jfc(r/3) is the Bessel function of the k th order and = 
eaEo. To obtain the current density, we expand £ s (p z )/^/o in Fourier series 
with coefficients £ rs 



£ s (pz, sAp ¥ 



Ss(pz)=J2 £ ™ eibrP *- ( 10 ) 



70 r/0 



Where 

a 



^ = o / £s(Pz)e- lbrp *dp z (11) 
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and expressing the velocity as 



v z (p z ,sA, Pip ) = = 70 Y,iar£ r .e ibr *', (12) 

dp z / — ' 



r^O 

the surface current density is then defined as 

2e 



(2tt/i) 2 



f{p)v z {p)dp. 
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or 



26 \ — f a 

jz = ( 2lT h) 2 ^ J ^ Pzi sAPv " ^w)^^*' sAp^dpz. (13) 

The integration is taken over the first Brillouin zone. Substituting Eqs. ([7]), 
Q and (10) into (13) we determined the current density for the zigzag CNT 
after averaging over a period of time as 

. _ 8e 7 o y- r y- f c V ^(rffl(l-i(nr + faj)T 
'■-^ata^^s'^S*^-^ l + «nr + fa,)r)» <14) 

When the electric field amplitudes are small (3 « 1, we can use the Bessel 
function approximation j± ~ (r/3/2) 2 and jo ~ 1 — r 4 /3 2 /4. Hence from Eqn. 
(14) we obtain real part of the current density to be, 

oo oo 

jz = jo ^ ^ T ^ ^ frs^rs 
r=l s=l 

(1 + (ut) 2 - (oj B rr) 2 ) 



[1 + ((w B r + w)r) 2 ][l + ((w B r + w)r) 5 



(15) 



and the imaginary part 



jz — jo ^ r ^ frs^rs 
r=l s=l 

x ( ^r(l - 3(rftr) 2 + (^r) 2 ) \ 

Ul + [lobttY] [1 + (( WjJ r + w)r) 2 ] [1 + (( WB r + ^)r) 2 ] J 1 j 

These are similar to the equations obtained in ref [25J for (3 « 1 . On 



the other hand we can re- write Eqn. ( 14 ) in the form of ref [26] as; 

oo oo n 

Jz = Jo £V ]T J 2 k {rp)j dc (E r + kE*)Y,frs£rs (17) 

r=l fc=— oo s=l 

Where E* = huj/a and for / 0, jdc(£lr + /ci?*) is the photon replicas. 
Expression ( JIT] ) is the superposition of the dc current density and its various 
photon replicas. We obtained the dc differential conductivity as; 

a z = (d Jz )/( d E o) = Jo ±r ± Ufa?) f;$rtkE*)$ t ^ 

r=l k=-oo V V ' > s =l 

(18) 

Where f2 = eaEo, (3 = eaE/ui for zigzag CNT and Q, = eaEo/^/3, (3 = 
eaEjypiuj for armchair CNT. This is exactly what is obtained by ref [27] 
except that the ID CNTs have the stark component (summation over r) 
which is absent in ref [27] and consequently making the ID CNTs more 
nonlinear than the 2D SL 
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3 Results, Discussion and Conclusion 

We present the results of a semiclassical theory of electron transport in a 
CNTs simultaneously exposed to both constant (dc) and ac electric fields. 
The Boltzmann's equation is solved in the frame work of momentum-independent 
relaxation time. An analytical expression for the current density is obtained 
for the situation /3 « 1. The nonlinearity is analyzed basically on the de- 
pendence of the current density on the frequency and the d.c electric field. 
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Figure 1: j z /jo — u/ojb curves for SL (thick solid curve) a-SWCNT (dotted 
curve) and z-SWCNT using r = 3 x l(T 12 s, T = 287.5if 

In Fig. [TJ we show the plot of the real part of the normalized current 
density (j z /jo) as a function of dimensionless frequency (co/ujb) for a super- 
lattice (thick solid curve), a a-SWCNT and a z-SWCNT, respectively. In all 
three systems, we observed that the real part of the differential conductivity 
is initially negative at zero frequency and becomes more negative with in- 
creasing frequency, until it researches a resonance minimum at a frequency 
just below the Bloch frequency (ojb t = 10) and then turning positive (res- 
onance enhancement) just below the Bloch frequency (see Fig. [I]). This 
resonance enhancement is indicative for terahertz gain at frequency without 
the formation of current instabilities induced by negative dc conductivity. 
It is worth to note that the effect is about twice and trice stronger in the 
a-SWCNT and z-SWCNT, respectively, in comparison to the SL, indicating 
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Figure 2: j z /j — E/E* curves for SL (thick solid curve) a-SWCNT (dotted 
curve) and z-SWCNT using r = 3 x l(T 12 s, T = 287.5A" 



stronger effective suppression of the current instability in CNTs. 

The dependence of j z on E a Eqn. (17) as well as a z Eqn. (18) on 



E a arising under the action of strong ac field. If ac field is applied along 
the CNTs axis in addition to the dc field, , new transport channels are 
opened at large ac field which are seen as distinctive peaks in the current 
density and the conductivity characteristics as shown in Figs. [2] and [3| 
respectively. We note a very steep positive slope in the neighborhood of the 
maximum electric field, which indicates that the differential conductivity is 
positive (PDC) in this region (see Fig. [2]). The PDC is considered one of 
the conditions for electric stability in the nonlinear system. The negative 
ac conductivity at the drive frequency appears when electric is weak. We 
observed that in the nonlinear regime with high-enough frequency, the dc 
differential conductivity is positive, while the large-signal high frequency 
differential conductivity remains negative (see Fig. [3]), which can be used 
for THz gain. 

The magnitude of the current density as well as dc differential conductiv- 
ity is an order of magnitude higher a-SWCNT and z-SWCNT in comparison 
to that of SL. The strong effects in CNTs are due to the higher density of 
free electrons in the CNTs. The mechanism of the nonlinearity in CNTs is 
due to the presence of the high stark components (summation with respect 
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to r, Eqns. (16 - 18) [|f5| which are absent in superlattice [2?] . 

In conclusion, we have used the semi-classical Boltzmann transport equa- 
tion to obtain the electric current density and the dc conductivity in a achiral 
CNT under the influence of simultaneous dc, ac electric field applied parallel 
to the tube axis. Our analysis on the current density versus electric field 
characteristics demonstrate negative high frequency differential conductiv- 
ity resonance enhancement, photon assisted peaks, suppression of current 
instability induced by negative dc conductivity and therefore the possibility 
of terahertz gain at room temperature without electric instabilities. 
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